Explanation of 100-fold reduction of spectral shifts for hydrogen on helium films 
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We show that helium film-mediated hydrogen-hydrogen interactions account for a two orders of 
magnitude discrepancy between previous theory and recent experiments on cold collision shifts in 
spin-polarized hydrogen adsorbed on a helium film. These attractive interactions also explain the 
anomalous dependence of the cold collision frequency shifts on the '^He covering of the film. Our 
findings suggest that the gas will become mechanically unstable before reaching the Kosterlitz- 
Thouless transition unless the experiment is performed in a drastically different regime, for example 
with a much different helium film geometry. 

PACS numbers: 67.63.Gh, 32.30.Bv, 67.25. bh, 32.70.Jz 



Two-dimensional (2d) spin polarized hydrogen is an 
exciting quantum fluid. It can support a superfluid- 
normal Kosterlitz-Thouless transition, which crosses over 
to Bose-Einstein condensation as the third-dimensional 
trapping is varied. Hydrogen experiments are unique in 
involving techniques from both low-temperature physics 
and atomic physics. The large particle numbers and light 
mass give a high Bose-Einstein condensation tempera- 
ture, and the interaction properties are known to high 
accuracy. 

In this context, it is remarkable that experiments on 
spin-polarized hydrogen adsorbed on a helium film dis- 
play a cold collision frequency shift which is two orders 
of magnitude smaller than theory predicts [l|. Here we 
show that a helium mediated interaction can explain this 
discrepancy and a similarly unexpected dependence of 
the shift on the concentration of '^He in the film. 

Since these experiments are in the "cold collision" 
regime 0], where the thermal de Broglie wavelength is 
longer than the interaction's effective range, they see ex- 
tremely sharp spectral lines which are shifted by inter- 
actions. These shifts are important: they limit hydrogen 
maser [1, Q and atomic fountain Q performance; the IS- 
2S spectral shift first demonstrated Bose-Einstein con- 
densation in hydrogen the Mott insulator-superfluid 
phase transition was probed in atomic gases in op- 
tical lattices 0, and the shifts revealed features 
of two-component Fermi gases across a Feshbach reso- 
nance Such spectra are promising probes of future 
strongly-correlated atomic systems. 

Calculating spectra. — Letting \ij) be the state with 
electron spin i and nuclear spin j, the relevant hydro- 
gen spin states in the experimental B = 4.6T field are 
|1) =cos6'|iT)-sin6'|Ti) and |2) = cos6i |n> +sin6l |iT) 
with negligible. Initially all of the atoms are in state 
Neglecting, for now, the helium film's degrees of 
freedom, the quasi-2d hydrogen's Hamiltonian is 
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mentum k and internal state j; e^^fc — k'^ /2m + Sj — 
the free dispersion of the effectively 2d gas, including the 
internal energy Sj and chemical potential fij] and Vij^k 
the Fourier space interaction potential between atoms in 
states i and j. Throughout we set h—1. 

The spectrum is measured by counting the number of 
atoms transferred from state |1) to state |2) after apply- 
ing a frequency u electromagnetic probe for a short time. 
In Ref. [ij , a combination of electron spin resonance and 
nuclear magnetic resonance drives the transition. In the 
rotating wave approximation the probe beam's Hamilto- 
nian [nl is i?p = npJ2k e"'^""'''""''^^^*V'2,fcV'i,fc + h.c. 

The radio frequency and microwave photons transfer neg- 
ligible momentum. 

Given that the range of the the potential (~ lA) is 
significantly less than the 2d interparticle separation (~ 
lOOA) the average absorbed photon energy shift is given 



in terms of the scattering amplitude [12 
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with 32(0 = (^I(r)^I(0)^i(0)^i(r))/(^t(0)V'i(0)) 
fij the i-j scattering amplitude evaluated at character- 
istic momenta of the system (here, set by the thermal 

de Broglie momentum), and cti = (tl)\{0)'ipi{0)^ the 2d 

density. In a dilute thermal gas 32(0) = 2. A similar 
formula holds in 3d, where the scattering amplitude is 
fij = A-Kaij with aij the 3d scattering length. In the case 
we consider, the spectrum is a delta function 0. 

In quasi-2d, where kinematics are 2d but the 3d scat- 
tering length is much less than the perpendicular con- 
finement length, we can construct a 2d interaction giv- 
ing identical low-energy scattering properties. Assum- 
ing harmonic confinement with oscillation frequency Wosc 



the effective 2d scattering amplitude is IJ, |l5|, ll6| 
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with T/jj ^ bosonic creation operators for states of mo- 



where qx = ^/rnkBT is the thermal momentum and I = 
^Jh/ [mujosc) is the length scale of z-axis confinement of 
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the hydrogen gas. Following Ref. [ij's discussion, the 
characteristic length scale for the confinement in their 
experiments is = — l/\/2mEa ^ 5 A where 

is the adsorption energy of the hydrogen on the helium 
film. The spectral shift is then 



5lj = 
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lo/ai2 - log{qp'^/2) 
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li a <^Iq this reduces to the simpler form 
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52(0)^X1. (4) 
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as for a 3d gas with density cti/Zq- 

The spectral shift in Eq. (O is one hundred times larger 
than experiment. We will show that this discrepancy is 
consistent with the corrections found by including the 
hydrogen-hydrogen interaction mediated by the helium 
film, Vincd- Although T4iod is state- independent, it alters 
5lj by reducing the probability that two particles are close 
enough to feel the spin dependent Vij . 

Helium film-mediated interaction. — Reference [l7| de- 
rived a hydrogen-only effective action assuming that (i) 
the helium film's excitations are non-interacting, (ii) 
the hydrogen-helium interaction falls off as l/r^, and 
(iii) that the hydrogen confinement length Iq is signif- 
icantly smaller than the hydrogen-helium separation 
Defining the total hydrogen density operator pt{p) = 
'^j'>pjip)4'j{p), the Fourier transform of Ref. [l3]'s ef- 
fective Hamiltonian is 

H,s = H-^IJ d'pdy pt(p)pt(p')K.cd(p - p') (6) 



with mediated pair interaction given by 



in terms of 



^0 = 
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with Aq defined as the parameter controlling the 
hydrogen-helium van der Waals interaction strength such 
that the potential is Vh-Ho = — (6Ao/n7r)r~^ with n 
the helium density, and with S = (6Ao(/)g/ri("*), A = 

\J MC'^ (217 + '^)' C distance from the hydrogen 
gas to the helium film, C3 the film's third sound speed, 
M the helium mass, c?o the film thickness, (3 the film sur- 
face tension, and 4>g = \/ndo (estimates of parameters in 
experiments are given later). The non-dimensionalized 
potential is defined as 



d^pd^p' A{p)A{p')Gt 
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where we have defined 
1 
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and GF{x) = Koix), (10) 
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with Kq the zero'th modified Bessel function. Note that 
A{p) = — (12(/)gAo/7rnC^)yi(p/C) is the hydrogen- helium 
van der Waals interaction at a separation p and Gf{p) — 
{n / 2dQTr I3)G F {p / is the helium film's Green's function. 

Eq. dS]) neglects Vmcd's frequency dependence, which 
arises from retardation: one hydrogen atom excites a he- 
lium excitation which propagates to another hydrogen 
atom. In a future work we will explore retardation ef- 
fects, but here we neglect them. 

We numerically compute V(^{x) as a function of x and 



^. Typical results are shown in Fig. 1(a) 



Spectral line shifts with mediated potential. — To evalu- 
ate spectral shifts via Eq. ([2]), we calculate the scattering 
amplitude of Vtot = Vij + Vmod- Since the bare hydrogen- 
hydrogen potential is short-ranged compared to Vmcd and 
the depth of Vmed will be sufficiently low, we may replace 
Vtot with V^^^ = V-^ + T4iod, where V-j is an arbitrary 
short range potential reproducing V^^'s scattering am- 
plitude. We numerically solve the resulting two-particle 
Schrodinger equation with potential T4ot, replacing the 
potential Vij with a boundary condition at small inter- 
particle separations. 

We use the following estimates in our calculations, 
taken from Ref.'s [i|,[ll[23: ~ 5l, at ^ .721, a^ = 
O.I7I, A ^ 50l, C3 - Im/s qr = ^Jmk BT ^ (30l)-\ 
8 ~ 0.265VMC|/2m, and C ~ SA. Fig.'s 



1(b) and 1(c) 



show the shifts for various parameter values, the princi- 
ple results of this paper. For visualization, we rounded 
out the resonances, where the dilute gas assumptions fail, 
by replacing f ^ [1/ f -\- ia\ ^ with a ~ 1. 



Fig. 1(b) shows the scattering amplitudes for the 1- 



1 and 1-2 scattering as a function of the mediated in- 
teraction strength, with all other parameters fixed at 
their typical values. The two vertical lines correspond 
to zero mediated potential and typical mediated poten- 
tial strengths. Fig. 1(c) shows the spectral shift as a 



function of the potential depth and the characteristic 
length A. The red circle indicates typical parameter val- 
ues, and the interior of the red box factor of 2 varia- 
tions of these. Fig. 1(c) shows that for typical values. 



the shifts are reduced by a factor of 7; simultaneously 
increasing potential depth and A by a factor of 2 larger 
than typical values (keeping other parameters typical) 
gives a 30-fold decrease. Similar reductions are possible 
when other "typical" parameters are individually varied 
(while holding all others constant): 2-fold variations in 
many of them can produce additional 2-fold reductions. 
Keeping in mind that the parameters themselves each 
could believably vary by factors of up to perhaps five, 
our calculation demonstrates that the helium-mediated 
interaction is a consistent and plausible source for the 
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frequency shift reduction. The agreement with experi- 
ment is even better when one considers that a 3d control 
experiment finds an ill-understood factor of 7 discrepancy 
in the 3d gas [2^ . Finally, one should note that ~ 40% of 
the change comes from simply using the more accurate 
formula of Eq. ^ in place of the approximate Eq. ([5]). 

^He film. — Our theory also accounts for unexpected 
effects of adding "^He to the film. It is known that adding 
^He reduces the hydrogen surface adsorption energy so 
one would naively expect an increase in the confinement 
length, decrease in density, and decrease in spectral shift. 
Instead, adding ^He to the superfluid helium film leads 
to an observed increase in the shift. 

Upon incorporating our theory of the helium surface, 
changing the adsorption energy Ea has two effects: it 
increases the confinement length / by ~2 — 10% and in- 
creases ^, the distance of the gas to the helium surface. 
We will assume that the fractional change in C is com- 
parable to the fractional change in I. Then Eq. ([8]) re- 
veals that the mediated potential Vmcd is proportional 
to cx l/l^, coming from 5, plus significantly weaker 
dependence from V'j's dependence on A through ^ — X/(. 

Thus, a 2 — 10% increase of I yields a '--^ 8 — 30% de- 
crease of the mediated potential strength. That such a 
decrease in mediated potential strength yields a signif- 
icantly larger spectral shift is clear from Fig. 1(b) the 



■^He doping shifts the typical parameters from the left 
dashed line to the right by 8 — 30%. The l/l factor in 
the spectral shift, due to the decreased 3d density, is in 
general insufficient to compensate. Using typical values 
for parameters yields a range of a 2-20% increase of fre- 
quency shifts upon adding '^He. This is consistent with 
the experimentally observed 25% increase, further evi- 
dence that the helium surface is playing a crucial role in 
the experiments. 

Validity of Approximations. — We have made a number 
of approximations. Here we enumerate the most impor- 
tant ones, and discuss which of them need to be addressed 
in future works through a more sophisticated theory. 

The most severe approximation we made is to follow 
Ref. [l7|, and neglect retardation in the induced poten- 
tial. Such effects are relevant when the phase velocity of 
the hydrogen excitations uj/k becomes large compared 
to C3 at characteristic velocities ^JksT /m and ener- 
gies T [2^. For r ^ 50 mK as in the experiments, 
w/A; ~ 5 m/s. In comparison, typical third sound veloci- 
ties are C3 ~ 1 m/s, so we expect retardation corrections 
may be significant. Including the frequency dependence 
in calculating the spectral shift is challenging, requiring 
solution of a coupled set of 2d partial differential equa- 
tions for each w. We expect that a sharp spectral peak 
survives, but with reduced spectral weight. 

Additionally, given the long length-scale of the me- 
diated potential, we should also critically examine the 
assumption that only the long-wavelength limit of the 
s-wave phase shift is needed to evaluate the cold col- 
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FIG. 1: (a) Rescaled mediated potential as a function, as 
given by Eq. ((9)1, with ^ = 1,4, 15 from top to bottom, (b) 
Scattering amplitudes / as a function of mediated potential 
depth Vd = (252)/(7r^A^MC|) for the triplet (dashed) and 
singlet (solid) scattering channels. Vertical lines indicate (i) 
the strength of the mediated potential using the typical pa- 
rameters given in the text and (ii) the absence of a mediated 
potential. The horizontal axis's units are converted assuming 
A = 50A, with A the helium film Green function's charac- 
teristic length, (c) Frequency shift reduction as a function 
of potential depth scale factor Vd = {2S^) /{n''^ XoMCi), with 
Ao = 50A and characteristic length A. Typical values of pa- 
rameters are shown by the red dot, while factor of two varia- 
tions comprise the interior of the red rectangle. Contour lines 
correspond to reductions of factors of 1/5, 1/15, and 1/30. 



lision frequency shift. Eq. requires that the areal 

— 1/2 

interparticle distance {n2£, ) is larger than the effec- 
tive range of the potential Re- Since Re ^ 250 A and 

— 1/2 ° 

n2£) ~ 100 — 300A, the approximation with only the 
fc = s-wave scattering shift may be sufficiently accurate. 
Similarly, the thermal wavelength is At ~ 30A while the 
cold collision regime strictly requires At ^ Re- This 
can smear out spectral lines somewhat, although since 
At ^ Re and the potential is rather shallow, the spectral 
lines may remain quite sharp. Here retardation helps us, 
as the slower moving atoms are the dominant contribu- 
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tors to the spectral peak, and these have a substantially 
larger de Broglie wavelength - a factor of 5 for the fastest 
contributing atoms [27 1 

Finally, we have neglected mass renormalization com- 
ing from virtual desorption- adsorption processes. This 
approximation is well-justified since the desorption pro- 
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cess is relatively slow 3, implying that the mass is 
renormalized by at most a few percent . 

Conclusions. — In summary, we have shown that incor- 
porating the helium surface into the theory of the hy- 
drogen gas can resolve the ~ 100-fold theoretical over- 
estimate of the spectral shifts compared to experiment, 
for a large parameter range consistent with experimen- 
tal estimates of the parameters. In addition, we have 
shown that adding ^He to the helium film increases spec- 
tral shifts, in agreement with experiment and in contrast 
to the naive theory. 

We have found that the helium surface induces an at- 
tractive interaction between the hydrogen atoms. It is 
useful to speculate on what other physical effects this 
interaction can cause. For example, can it drive a me- 
chanical instability? While thermal pressure will stabilize 
the gas at the temperatures currently being studied, we 
believe that as the temperature is lowered that a "col- 
lapse" might occur, similar to the ones seen in atomic 
gases 21,22, [23| • Importantly, the Kosterlitz-Thouless 
transition requires repulsive interactions, so the medi- 
ated interaction may make eliminate the possibility of 
achieving superfluidity without significantly altering ex- 
perimental parameters. 

Our calculation suggests an experiment which would 
further constrain the modeling of this system. Replacing 
hydrogen with deuterium or tritium drastically changes 
the parameters in our theory by increasing the binding 
energy to the helium film, by factors of roughly two and 
three, respectively [3]. However, the form of the theory 
remains the same, so one can compare the expectations 
on the basis of our theory with experiment. 

Finally, we point out possible ramifications of our the- 
ory to ultracold atomic gases and elsewhere. The key to 
our findings is that in quasi-2d, infinitesimal attractive 
interactions generate scattering resonances. These are 
pushed to finite depth by the 3d interaction, but a gen- 
uine 2d interaction easily overwhelms the effects of the 3d 
interaction. One can imagine our present theory apply- 
ing to quasi-2d fermi-bose or bose-bose mixtures where 
one system mediates an interaction for the other species. 
Interestingly, similar physics might arise in layered sys- 
tems: one layer effectively mediates an interaction in ad- 
jacent layers - incorporating these effects in a consistent 
manner could lead to dramatically modified interaction 
properties for each layer. 
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